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We discuss the time evolution of the wave function which is solution of a stochastic 
Schrodinger equation describing the dynamics of a free quantum particle subject to sponta- 
neous localizations in space. We prove global existence and uniqueness of solutions. Observ- 
ing that there exist three time regimes, namely the collapse regime, after which the wave 
function is localized in space; the classical regime, during which the collapsed wave function 
moves along a classical path and the diffusive regime, in which diffusion overlaps significantly 
the deterministic motion we study the long time behavior of the wave function. We assert 
that the general solution converges a.s. to a diffusing Gaussian wave function having a finite 
spread both in position as well as in momentum. This paper corrects and completes earlier 
works on this. 



INTRODUCTION 



Stochastic differential equations (SDEs) in infinite dimensional spaces are a subject of growing 
interest within the mathematical physics and physics communities working in quantum mechanics; 



they are currently used in models of spontaneous wave function c ollap se QsjO, J, 5,^,13, 8, |^ 
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in the theory of continuous quantum measurement [l5l . 



and in the theory of open quantum systems 
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26( 1 ■ In the first case, the Schrodinger equation 



is modified by adding appropriate non-linear and stochastic terms which induce the (random) 
collapse of the wave function in space; in this way, one achieves the goal of a unified description 
of microscopic quantum phenomena and macroscopic classical ones, avoiding the occurrence of 
macroscopic quantum superpositions. In the second case, using the projection postulate, stochastic 
terms in the Schrodinger equation are used to describe the effect of a continuous measurement. In 
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the third case, shghtly generahsing the notion of continuous measurement to generic interactions 
with environments, SDEs are used as phenomenological equations describing the interaction of a 
quantum system with an environment, the stochastic terms encoding the effect of the environment 
on the system. Looking directly at the stochastic differential equation for the wave function, 
rather than the deterministic equation of the Lindblad type for the statistical operator has some 
advantages with respect to the standard master equation approach, e.g. for faster numerical 
simulations [2^. 

Among the different SDEs which have been considered so far, the followin g eq uation, defined 



in the Hilbert space Ti = C' 

n2 



d^Jt 



, is of particular interest 
A 



16|, 



29 



3Q|, 



31 



32 



33, 



34 



-l^'^t + ^{q-{q)t)dWt - ^{q-{q)t?dt 



351, m,\m- 

i'o = i'- (1) 



The first term on the right-hand-side represents the usual quantum Hamiltonian of a free particle 
in one dimension, p being the momentum operator. The second and third terms of the equation, 
as we shall see, induce the localization of the wave function in space; q is the position operator 
and {q)t denotes the quantum expectation {tptlq'ipt) of q with respect to ^pt■ The parameter A is a 
fixed positive constant which sets the strength of the collapse mechanism, while Wt is a standard 
Wiener process defined on a probability space {^},J^,¥) with filtration {J-t,t ^ 0}. 

Eq. ([1]) plays a special role among the SDEs in Hilbert spaces because it is the simplest exactly 
solvable equation describing the time evolution of a non-trivial physical system. Within the theory 
of continuous quantum measurement, it describes a measurement-like process designed to measure 
the position of a free quantum particle; within decoherence theor y it represents one of the possible 
unravellings of the master equation first derived by Joos and Zeh [38t | . Within collapse models (like 



GRW-models), it may describe the evolution of a free quantum particle (or the center of mass of 
an isolated system) subject to spontaneous localizations in space 2] in the following sense. 
Realistic models of spontaneous wave function collapse are based on a more complicated stochastic 
differential equation: The difference between Eq. (fTl) and the equations of the standard localization 

n n 

models such as GRW [1] and CSL [2] is most easily described on the level of the Lindblad equations 

t){tpt\], induced by the stochastic dynamics of 



(2) 



for the respective statistical operators pt := 
the wave function. By virtue of Eq. ([T|) (see e.g. 



d 

dtp' 



« r 2 ^ A 



with the "Lindblad term" in position representation 



(x - yfpt{x,y) . 



(3) 



For the GRW dynamics as described in [l| the corresponding Lindblad term of the GRW master 
equation in the position representation reads: 



Ar 



1 



-a(a;-s/)V4 



Pt{x,y) . 



(4) 



When the distances involved are smaher than the length Xj \/a ~ 10 ^ cm characterizing the model 
we have that 



Ac 



1 



-a(x--j/)2/4 



A 



-[X - y) 



for: \x — y\<^ Xj ^fa . 



(5) 



Accordingly, the stochastic dynamics of Eq. ([T]) approximates-at least on the statistical level-the 
GRW dynamics for all atomic and subatomic distances. Since this is a regime of growing interest 



4QI, 



41 



42l ] it is reasonable to study now first the simpler equation Eq. ([T]). 



Eq. ([T|) is non-linear. Non-linearity is a fundamental ingredient because only in this way it 
is possible to reproduce the collapse of the wave function. It is well known how to "linearize" 
the equation, i.e. how to express its solutions as a function of the solutions of a suitable linear 



SDE 



29|, 



43( 1 . We briefly review this procedure. 



Let us consider the following linear SDE: 



dt + VXqdCt - ^q^dt 
n 2m 2 



(6) 



defined in the same Hilbert space TC = £^(M); the stochastic process is a standard Wiener 
process with respect to the probability space (fi,.?-", Q) and filtration {J^tjt > 0}, where Q is a 
new probability measure whose relation with P will soon be established. This equation does not 
conserve the norm of the state vector, as the evolution is not unitary; we therefore introduce the 
normalized state vectors: 



(7) 



otherwise; 

A standard application of Ito calculus shows that, if (j)t solves Eq. ([6]), then ipt defined in ([Tj) solves 
the following non-linear SDE: 

.2 _ _ 



-{^dt + V\{q - {q)tmt - 2VX{q)tdt) - ^{q - {q)tfdt 



(8) 



for the same initial condition ip = (p. 

Eq. (I8|) is a well defined collapse equation, however it is not suitable for physical applications, as 
the collapse does not occur with the correct quantum probabilities. This can be seen by analyzing 
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the time evolution of particular solutions, such as Gaussian wave functions; it can also be easily 
understood by noting that there is no fundamental difference between Eq. ([8]) and Eq. Q, since 
any solution of Eq. ([8]) can be obtained from a solution of Eq. ([6]) simply by normalizing the wave 
function. In turn, Eq. ([6]) does not contain any information as to why the wave function should 
collapse according to the Born probability rule, i.e. the Wiener process is not forced to pick 
most likely those values necessary to reproduce quantum probabilities, during the collapse process. 

The way to include such a feature into the dynamical evolution of the wave function is to replace 
the measure Q with a new measure (which will turn out to be the measure P previously introduced) 
so that the process ^t, according to the new measure, is forced to take with higher probability the 
values which account for quantum probabilities. This is precisely the key idea behind the original 
GRW model of spontaneous wave function collapse [l| : the wave function is more likely to collapse 
where it is more appreciably different from zero. The mathematical structure of the GRW model 
suggests that the square modulus \\4>t\\'^ should be used as density for the change of measure. We 
now formalize these steps. 



In 



29|, Holevo has proven that \\4>t\\ is a martingale satisfying the equation: 



Jo 



(9) 



when 1 1 00 IP = li and from now on we will always assume that this is the case, WcptW^ can be used 
as a Radon-Nikodym derivative to generate a new probability measure P from Q, according to the 
usual formula: 

F[E] := EqilEUtf] \/ E £ J^t V t < +oo, (10) 

where l^: is the indicator function relative to the measurable subset E. The martingale property, 
together with the property EQ[||0i|p] = 1, guarantee consistency among different times, so that (jlOp 



defines a unique probability measure P. In the following, for simplicity we will write dF/dQ = \\4't\\'^- 
One can then show that Eq. dH), with the stochastic dynamics defined on the probability space 
(0,jr, P) in place of (r2,J^, Q), correctly describes the desired physical situations. 

A drawback of the change of measure is that the equation is defined in terms of the stochastic 
process .^t, which is not anymore a Wiener proves with respect to the measure P, as it was with 
respect to the measure Q. This can be a source of many difficulties, e.g. when analyzing the 
properties of the solutions of the equation. The disadvantage can be removed by resorting to 
Girsanov's theorem, which connects Wiener processes defined on the same measurable space, but 
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with respect to different probability measures. According to this theorem, the process 

Wt := it - 2VA f{q)sds, 
Jo 



(11) 



is a Wiener process with respect to {il.,J-,¥) and filtration {J-t,t > 0}, and thus is the natural 
process for describing the stochastic dynamics with respect to the measure P. It is immediate to 
see that, once written in terms of Wt, Eq. ([8]) reduces to Eq. ([1]), thus the link between Eq. ^ 
and ([ID is established. The above discussion should also have given a first idea of why SDEs like 
Eq. (dl) are those which are used in Quantum Mechanics to described the collapse of the wave 
function; we will come back on this point later in the paper. 

The first important problem to address concerns the status of the solutions of Eq. 



In 



291, 



Holevo has proven the existence and uniqueness of topological weak solutions of a rather general 
class of SDEs with unbounded operators, to which Eq. ^ belongs. (See the end of the section 
for the notation.) The problem of the existence and uniqueness of topological strong solutions 



of Eq. dll) has been addressed in 



281]; there however, the proof relies on the expansion of wave 



functions in terms of Gaussian states, which in general is problematic and requires special care, as 
shown in [44 1 . An explicit representation of the strong strong solution of Eq. ^ has been given 



m 



351 ]: the representation is written in terms of path integrals and is not particularly suitable for 



analyzing the time evolution of the general solution. A much more convenient representation, given 
in terms of the Green's function of Eq. ([6]), has been first derived in 
reads: 



331 ]: the Green's function 



at 



— {X +y )+ fitxy + atx + hy + ct 



Gt{x,y) = Ktexp 
the coefficients Kt, at and Pt are deterministic and equal to 

Kt 
at 



A 



VTT sinh vt 

2A , 

— cothut, 

V 

I3t = 2 - smhr^vt, 

V 



while the remaining coefficients are functions of the Wiener process ^4: 

at = ^fXsnihr^vt / sinh us , 
Jo 



^ _ 2ihX 
* m V Jq 
ih 



sinh vs 



ds, 



Ct 



m Jo 



at ds. 



(12) 

(13) 
(14) 
(15) 

(16) 
(17) 
(18) 
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In the above expressions, we have introduced the following two constants: 

„ . . . 2^-. (19) 

As we shall see, the parameter uj, which has the dimensions of a frequency, will set the time scales 
for the collapse of the wave function. The representation in terms of the Green's function ()12p . as 
we said, is particularly suitable for analyzing the time evolution of the general solution of Eq. ([B]), 
and thus of Eq. ([T|), even though we will see that, when studying the long time behavior, another 
representation is more convenient. 

Our first result concerns the meaning of the solution of Eq. ([6]) in terms of 

Mx) ■■= J dyGt{x,y)4>{y) (20) 

for given initial condition 0. 

Theorem 1 (Solution): let 4)t be defined as in (pOj) : then the following three statements hold 
true with Q-probability 1: 

1. (/> E C'^{W) (Pt e C\R), (21) 

2. (p £ ^b(^) =^ is a topological strong solution of Eq. (22) 

3. ^ e C'^iR) lini||(/)t-(/)|| =0, (23) 

where £^ (M) is the subspace of all bounded functions of (M) . 

Having the explicit solution of the Eq. ([6]), and thus of Eq. ([T|), the next relevant problem is 



14 



36( 1 have shown 



to unfold its physical content. Previous analysis of similar equations 
that one can identify three regimes, which are more or less well separated depending on the value 
of the parameters A and m. 

1. Collapse regime: A wave function having an initial large spread, localizes in space, the 
localization occurring in agreement with the Born probability rule. 

2. Classical regime: The localized wave function moves in space like a classical free particle, 
since the fluctuations due to the Wiener process can be safely ignored. 

3. Diffusive regime: Eventually, the random fluctuations become dominant and the wave func- 
tion starts to diffuse appreciably. 
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It is not an easy task to spell out rigorously these regimes and their properties. We shall however 
be a bit more specific on this in the following section. We shall afterwards focus on the simplest 
regime, namely the diffusive one, which in fact has been intensively looked at in the previous 
years [7j, 



17 



24 



32 



33 



36| and we shall prove a remarkable property of the solutions of Eq. ([T]): 



Any solution converges almost surely to a Gaussian state wave function having a fixed spread. 

Theorem 2 (Large time behavior): let ipt be a solution of Eq. then under conditions 
which we will specify, the following property holds true with P-probability 1: 



lim ||-0t 
t— >oo 



0, 



(24) 



where defined in (jll6p . is a Gaussian wave function with a fixed spread both in position and 
momentum. 



Theorem 1 and Theorem 2 have been extensively discussed before in the literature 



33, 



341, 



36l |. proving that the community has devoted much attention to the problem. However, 



these proofs are not complete or flawed. Concerning Theorem 1, in particular Statement 3 was not 



proven 



28|, 



32 



33, 



34l |. While Statements 1 and 2 are rather straightforward conclusions from the 



Gaussian kernel of the propagator, the third Statement is much more subtle and does not follow 
from purely analytical arguments. Concerning Theorem 2, none of the previous proofs is decisive. 
In Q, Q> 

the major flaw was that it was overlooked that the eigenfunction expansion of the 



relevant dissipative operator (not self-adjoint) does not give rise to an orthonormal basis. In [l7l ]. 
the long time behavior was analyzed by expanding the general solution in terms of coherent states, 



while in [24 
solution; in 



it was analyzed by scrutinizing the time evolution of the spread in position of the 
44!] it has been shown that both approaches are not conclusive. Finally, [7] proposed 



Theorem 2 as a conjecture, but shows stability of "0*° only against small 



perturbations. Building 



45l | and 46| the general 



on previous work of Holevo, Mora and Rebolledo recently enhanced in 
theory of stochastic Schrodinger equations. In particular they developed criteria for the existence 
of regular invariant measures for a large class of stochastic Schrodinger equations as an important 
step towards an understanding of the large time behavior. Until now however the only complete 
and detailed results on the large time behavior seem to be Theorems 1 and 2. 

We conclude this introductory section by summarizing the content of the paper. In Sec. HI] we 
will present a qualitative analysis of the time evolution of the general solution of Eq. ([T]); we will 
discuss the three regimes previously introduced, giving also numerical estimates, and we will set 
the main problems which we aim at solving. In Sec. lIIII we will analyze the structure of the Green's 



8 



function (jl2p and prove theorem 1. In Sec. IIVI we will introduce another representation of the 
general solution of Eq. ([T]), which is more suitable for analyzing its long time behavior. Sec. IVl will 
be devoted to the proof of theorem 2. Finally, Sec. IVII will contain some concluding remarks and 
an outlook. 

Notation. We will work in the complex and separable Hilbert space vC^(M), with the norm and 
the scalar product given, respectively, by || • || and (•!•). We will also consider the subspace 'C^(M) 
of all bounded functions of £^(M). Given an operator O, we denote with I'(O) its domain and with 
71(0) its range. 

Since in some expression the real and imaginary parts of some coefficients appear, we introduce 
for ease of readability the symbols or will denote the real part of the complex number z, 
while z^ or Zi will denote its imaginary part. 

Given the linear SDE ([6]) , a topological strong solution is an £^-values process such that for any 
t > 0, 

rt J2 rt \ rt 



[ ^(t)sds + \/A / q(t)sdWs - ^ [ q'^cpsds (25) 
Jo 2m Jo 2 Jo 



holds with Q-probability 1. A topological weak solution instead is an £^-values process such that 
for any t > and for any x ^ T^ip^) ^ T>{q^)^ 



{x\4>t) = (xl0) 



I f ^{pW(t>s)ds + VA f {qx\(t>s)dWs - \ l\qW<Ps)ds (26) 
n Jo 2m Jo 2 Jo 



holds with Q-probability 1. Topological strong and week solutions for the nonlinear SDE ([T]) are 
defined in a similar way. 

There is also a distinction between strong and weak solutions in a stochastic sense {4?! , depending 
on whether the probability space, the filtration and the Wiener process are given a priori (strong 
solution) or whether they can be constructed in such a way to solve the required SDE (weak 
solution). Throughout the paper we will deal only with strong solutions in the stochastic sense. 



II. TIME EVOLUTION OF THE GENERAL SOLUTION 



We begin our discussion with a qualitative analysis of the time evolution of the general solution 
of Eq. ([T]); we will spot out the regimes we introduced in the previous section, corresponding to 
three different behaviors of the wave function. These regimes of course depend on the value of the 
mass m of the particle and also on the value of the coupling constant A which sets the strength of 



the collapse mechanism. As discussed e.g. in [36], it is physically appropriate to take A proportional 
to the mass m according to the formula: 



A := An 



m 



mo 



(27) 



where Aq is now assumed to be a universal coupling constant, while rriQ is taken equal to the mass of 
a nucleon (~ 1.67 x 10^^^ kg). To be definite, in the following we take Aq — 1.00 x 10~^ m~^ sec~^, 
so that the localization mechanism has the same strength as that of the GRW model [1^. Though, 
as we discussed in the introduction, Eq. ([1]) is used also in the context of the theory of continuous 
measurement as well as in the theory of decoherence, for brevity and clarity in the following we 
will only make reference to its application within models of spontaneous wave function collapse. 

1. The collapse regime. The first important effect of the dynamics embodied in Eq. ^ is 
that a wave function, which initially is well spread out in space, becomes rapidly localized. This is 
most easily seen through the Green's function representation of the solution. The Green's function 
Gt{x,y) in (fT2]) can be rewritten as follows 



Gt{x,y) = Kt exp 
where we have introduced the new parameters: 



— X +atx + ct 



exp 



(28) 



at 



at 



Ct 



Yt" 



at 



at + 



ct + 



PI 
at 

f3tbt 
1? 



2A 



tanh vt, 



2at' 
Ptx + bt 

at 



(29) 
(30) 
(31) 
(32) 



The y-part of Gt{x,y) is a Gaussian function whose spread in position (equal to l/y^af) rapidly 
decreases in time, and afterwards remains very small. In particular, we have: 



r 2 



2A sinhwt — sintijt 



a^ 



cosh ojt — cos ujt 



- At ~ (3.99 X lO^^m-^ Kg"^ sec-i)mt t < uj- 
3 



2A 

— ~ (2.39 X lO^^m"^ Kg-i)m t +oo. 



(33) 



with UJ ~ 5.01 X 10~^ sec~^ independent of the mass of the particle. 

Let us introduce a length i, and let say that a wave function is localized when its spread is 
smaller than £. For sake of definiteness, we take i ~ 1.00 x 10^'' m, corresponding to the width of 
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the collapsing Gaussian of the GRW model. By means of this length, we can define the collapse 
time ti as the time when the spread of the y-part of the Green's function Gt{x, y) becomes smaller 
than i. By using the small time approximation of af given in ()33p . we can set: 

3 2.51 X IQ-" Kg sec 

As we see, and as we expect, this time decreases for increasing masses, i.e. for increasing values of 
A, and is very small for macroscopic particles. 

Let us assume that the initial state (p{x) is not already localized, and in particular that it does 
not change appreciably on the scale set by i; this is a physically reasonable assumption when (p 
represents the state of the center of mass of a macroscopic object. In this case, from the time ti 
on, the y-part of the Green's function Gt{x,y) acts like a Dirac-delta on and the solution at 

time t of the linear equation can be written as follows: 



[2n 

(t^tix) ^ \ — Ktexp 
V a* 



— x +atx + ct 



HYf); (35) 

This is a Gaussian state whose spread is controlled by at, which evolves in time in a way similar 
to at; in particular: 

' 2Xt ~ (1.20 X lO^^m"^ K:g-i sec-i)mt t < 
„ 2 X sinh ujt + sin ujt I . „. 

«f = . = { (36) 



UJ cosh Lot + cos LUt 



2A 

— ~ (2.39 X 1029m-2 Kg-^)m t +oo. 

UJ 



As we see, the spread is well below £, for any t >ti. We can the conclude that, for times 

greater than the collapse time, any state initially well spread out in space is mapped into a very 
well localized wave function. 

An important issue is where the wave function collapses to, given that the initial state is spread 
out in space. We now show that the position of the wave function after the collapse is distributed 
in very good agreement with the Born probability rule. 

A reasonable measure of where the wave function is, after it has collapsed, is given by the 
quantum average of the position operator {q)t. Accordingly, the probability for the collapsed wave 
function to lie within a Borel measurable set A of M can be simply defined to be P^°"[A] := ¥[uj : 
{q)t G A]. Though this probability is mathematically well defined for any Borel measurable subset 
A, it is physically meaningful only when A represents an interval A much larger than the spread 
of the wave function itself, or a sum of such intervals. In such discussed in [148], one can 

show that: 

Ff\A] ~ Ep[||PAV'tf] = ! Pt{x)dx, (37) 

J A 
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where Pa{x) is the characteristic function of the interval A of the real axis and pt = Ep[|'0t(x)p]. 
The idea behind the approximate equality ()37p is that when ipt lies within A, then PA'fpt — ipt, so 
that llPAV'tP is almost equal to 1, while when it lies outside A, it is practically 0. The critical 
situations, which require special care, are those when the wave function lies at the edges of A. 



In 



3a] it has been proven that: 



Ptix) 



^ I dye-^^y^n^ + y), = |^ (2.27 x lO^m^^ Kg-^ sec^) ^, (38) 



where pf'''^{x) = \'ipf'''^{x)\'^ and ip^'''^{x) is the solution of the standard free-particle Schodinger 
equation, for the given initial condition (p{x). For the times we are considering (t = ti), the 
Gaussian term in (j38p is much more peaked than any typical quantum probability distribution 
Pt'^'^ix), and consequently acts like a Dirac-delta on it; accordingly, pt{x) ~ pf''''(x). Finally, for 
macroscopic systems and for the times we are considering, the wave function solution of the free- 
particle Schrodinger equation does not change appreciably, implying that pf^^^x) ~ Pq'^^x) = 
|(/)(a;)p, which means precisely that the collapse probability is distributed in agreement with the 
Born probability rule. 

2. The classical regime. After time ti, we are left with a wave function which, when m is 
the mass of a macroscopic particle, is very well localized in space, almost point-like. This is the 
way in which collapse model reproduce the particle-like behavior of classical systems, within the 
framework of a wave- like dynamics. The relevant question now is to unfold the time evolution of 
the position and momentum of the wave function, to see whether it matches Newton's laws. 

When the wave function is well localized in space (t > ti), one can reasonably assume that it can 
be approximated with the Gaussian state to which — as we shall see — it asymptotically converges 
to. We will analyze the time evolution of such a Gaussian state in the following, and we will see 
that its mean position xt and momentum hkf evolve in time as follows (see Eqs. (|142p and (jl43p ): 

Xt = Xt, + -ktAt - ti) + VX- [ Wsds + J -{Wt - Wt,), (39) 
m ITT- Jti \ m 

kt = kt,+^{Wt-Wt,). (40) 

We can easily recognize in the deterministic parts of the above equations the free-particle equations 
of motions of classical mechanics describing a particle moving along a straight line with constant 
velocity; the remaining terms are the fluctuations around the classical motion, driven by the 
Brownian motion Wt- The important feature of the above equations is that these fluctuations, for 
macroscopic masses, are very small, for very long times. As a matter of fact, if we estimate the 
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Brownian motion fluctuations by setting Wt ~ we have for the stochastic terms in Eq. ([39]): 

VA- t Wsds ^ - VA- t3/2 ^ (1.63 X 10-22m Kg"/' sec^^/S) (41) 
m 3 m -^m 

.[^{Wt -WtA ~ a/^ ~ (1.02 X 10-1 V sec-i/2) (42) 

V m V rra V m 

We see that the random fluctuations decrease with the square root of the mass m of the particle, 
which means that the bigger the system, the more deterministic its motion. This is how collapse 
models recover classical determinism at the macroscopic level, from a fundamentally stochastic 
theory. 

We can introduce a time t2, defined as the time after which the fluctuations become larger than 
L; we can set e.g. L ~ 1.00 x 10^^ m. Since the fluctuations in ([^T|) grow faster as those in (jH]), 
we can set: 

/ 3 L Til \ 

t2 ~ ( ~ (3.55 X 10^2 sec m"i/^)^ ~ (1.13 x 10^ y m'^/^)^. (43) 



2^h 

The time ^2 defines the time interval [ti,t2] during which the classical regime holds. As we see, for 
macroscopic systems this is a very long time, much longer than the time during which a macro- 
object can be kept isolated from the rest of the universe, so that its dynamics is described by 
Eq. dH). 

To summarize, during the classical regime, which for macroscopic systems lasts very long, the 
wave function behaves, for all practical purposes, like a point moving deterministically in space 
according to Newton's laws. In other words, the wave function reproduces the motion of a classical 
particle. 

3. The diffusive regime. After time t2, two new effects become dominant: First, the wave 
function converges towards a Gaussian state, as we shall prove. Second, the motion becomes 
more and more erratic: the dynamics begins to depart from the classical one, showing its intrinsic 
stochastic nature. 

A thorough mathematical analysis of these time regimes and their main properties is still lacking. 
In this paper, as we have anticipated, we focus now only on the long time behavior of the solutions 
of Eq. ([T|), leaving the study of the remaining properties as open problems for future research. 



III. SOLUTION OF THE EQUATION 



In the first part of this section we derive the Green's function (jl2p in a way which will make 
clear the connection between Eq. ([6|) and the equation of the so called non-self-adjoint (NSA) 
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harmonic oscillator 5l|, |52| . l53l |. This connection is important for two reasons; from a physical 
point of view, it will bring a deep insight on how the collapse of the wave function actually works. 
From a mathematical point of view, it will allow to prove rigorously both the theorem 1 and 2 
presented in the introductory section. 

A way to connect Eq. Q with that of the NSA harmonic oscillator is to apply suitable transfor- 
mations to the wave function in such a way to transform the SDE in a Schrodinger-like equation. 
We will do this in two steps. We present this section in detail for convenience although the approach 
goes back to Kolokoltsov js^. 



1. Reduction of Eq. ([6]) to a linear differential equation with random coefficients. 
The idea is to remove the stochastic differential term \/\qd^t from Eq. ([6]): borrowing the language 
of quantum mechanics, we shift to a sort of interaction picture by defining a suitable operator 
which maps the solution of Eq. ^ to the solution of a new equation which does not have that 
stochastic term. To this end, let us consider the operator Qa : ^{Qa) C £^(R) defined 

as follows: 

Q„0(x) = e'^Xx), aeC; (44) 

where V{Qa) is defined as the set of all 4>{x) £ C'^iR) such that e^'^Xx) G C^{R). It should be 
noted that, in general, the operator Qa is unbounded and its domain V^Qa) is dense in £^(]R) but 
not coincide with it. We will settle al technical issues in the second part of the section. We now 
define the vector: 



an easy application of Ito calculus shows that satisfies the differential equation: 

o2 



dt, 



(45) 



(46) 



The stochastic differential ^/Xqd^t has disappeared; in turn, the free Hamiltonian /2m has been 
replaced by the operator Q_^^^ {p^ /2m) Q~^^ which, due to the specific commutation relations 
between q and p, takes the simple form: 



Q 



Eq. ()46p can then be re-written as follows: 



dt * 



P 

2m 



ihX ( 



ih 



m 



Xh^ 

2m 



(47) 



(48) 
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This is a standard differential equation with random coefficients; note that the operator on the 
right hand side is not self-adjoint, due to the presence of the second and third term. The last term 
of Eq. (l48]) is a multiple of the identity operator and can be removed by defining: 



we then obtain: 



dt * 



exp 



ihX I 



ih 



tp 



(49) 



(50) 



2m m 

The third term on the right-hand-side contains a time dependent coefficient, and the next step 
aims at removing it. 



2. Reduction of Eq. ([50|) to a differential equation with constant coefficients. The 
idea we now follow is to perform a transformation similar to a boost. We introduce the operator 
Va : V{Va) ^ /:^(M) ^ C'^{R) defined as: 



a G C, 



(51) 



where T>{Va) is the set of all (j){x) £ C'^(M) which can be analytically continued to the line x + a 
in the complex space C, and such that (j){x + a) E £^(M). Similarly to Qa, also Va is in general an 
unbounded operator and its domain I?('Pa), though being dense, does not coincide with £^(M); we 
will come back to this point later in this section. We define the operator: 



(52) 



where the coefficients Of, bt and Cf, yet to be determined, will turn out to be complex random 
functions of time. One can easily verify that: 



VtqV, 



-1 



+ h 



VtpV^ ^ = p + ct, 



and similarly for higher powers of q and p. Let us define the vector: 



(53) 
(54) 

(55) 



which solves the equation: 



7. ^ 



ihXq 

2m 



; 1 rrA 

h Q H VA^t p- 

m m I 



+ {ct-2ih\bt)q + [at + Ct bt + 



1 
2m 



ih 



m 



V\£,tCt - ihXbl 



(56) 
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The time-dependent part of the equation can be removed by requiring that at, bt and q satisfy the 
first-order differential equations: 



mbt — ct = —ih-\/~\^t bo = 0, 

ct-2ihXbt = Co = 



(57) 



and 



1 ih 

ht + ihXhi + —cl y/XitCt = 0, ao = 0. (58) 

2m m 

The first two equations form a non-homogeneous hnear system of first order differential equations, 
which has a unique Q-a.s. continuous random solution; the third equation instead determines the 
global factor a^, which is also random. With such a choice for the three parameters, Eq. (j56p 
becomes: 



7? 

— ihX q 

2m 



ipt, (po = 4>, (59) 



which is the equation of the so-called non-self-adjoint (NSA) harmonic oscillator, whose solution 
and most important properties are well known. Before continuing, we note that in the case of 
a more general Hamiltonian H = /2m + V{q) appearing in Eq. ([1]) in place of just the free 
evolution /2m, the potential V{q) would have been transformed, when going from Eq. ()50p to 
Eq. ([59j) . according to the rule: VtV {q)Vt^ = V{q+bt); in this case, we would not be able to remove 
completely the time-dependent terms from the equation and we would not be able to reduce the 
original equation to one, whose solution is known. However, besides the free particle case, all 
equations containing terms at most quadratic in q and p (among them, the important case of the 
harmonic oscillator) can be solved in a similar way. 

The solution of Eq. ()59p admit a representation in terms of the Green's function: 



GT\^,y) = \l ^exp 

" VTTsmhvt 



- — (x^+y'^)coth.vt + 2 — xysmh. ^ vt 

V V 



(60) 



with V and uj defined as in ()19p . In this way we have established the link between the solutions of 
the SDE (l6|) and those of the equation for the NSA harmonic oscillator (|59p . which we summarize 
in the following lemma, whose proof is straightforward. 

Lemma lllli l: Let T^^sa ^j^g evolution operator represented by the Green's function Gf^^{x,y) 
and % the one represented by Gt{x,y); then: 

% ^ e^p{m/h) Q^5,+(,,,/,)T'„,,,/,7;-^^ (61) 
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where the two random functions bt and q solve the hnear system (j57p . and ^t, which includes all 
global, i.e. independent of x, phase factors, solves the equation: 



2m m 2m 



= 0. 



(62) 



We now proceed to prove in which sense (pt '■= %:4> is the topological strong solution of Eq. ^ 
for the given initial condition (p. We first need to set some properties of the Green's function 
G^^^{x, y) which will be necessary for the subsequent theorem. 



Lemma [Tm 2: The absolute value of G^^^{x,y) is equal to: 



2A 



y vrwV cosh u>t — cos u>t 
where we have introduced the following quantities: 

sinhwt — sin 



exp 



-—[x +y )pt + 4: — xyqt 



UJ 



Pt 

Qt 



coshcjt — coswt ' 

sinh 01^/2 cos ujt/2 — cosh cot/ 2 sin ujt/2 



(63) 

(64) 
(65) 



cosh Ldt — COS tot 

note that the function pt is positive for any t > 0. The integral of |G^^*(x,y)p with respect to y 
is equal to: 



/ 



2X 



7ruj{sinh ujt — sinwt) 



exp 



^ Pt 



(66) 



A simple calculation shows that pi — Aq^ > for any t > 0; this means that Gf^*(x, •), taken as a 



function of y, belongs to £'^(M) for any x G M and t > 0; moreover: 



dx\\Gf^{x,-)f < +00 



for any t > 0. 



(67) 



Finally, the following expression holds true: 

dy\e'^Gr^{x + a,y)\^ = ' 



with 



2A r A 

— : ^ exp < — 2— 

TTU) [smh Lot — smut) [ uj 



Pj - ^2 

Pt 



+ 2 ptOR + Pi-tti - 4 h 26r 1 X 



Pt 



Pt 
Qt 



+ pt{al - of) + 2p^as,ai - 4 



sinhwt + sinwt 



Pt 



coshcjt — coswt ' 

sinh ujt/2 cos /2 + cosh cot/ 2 sin cjt/2 
cosh cot — cos a;t 



(68) 

(69) 
(70) 
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The above formulas imply that, for any a, 6 E C, for any x G M and for any t > 0, the function 
^bxQNSA(^^ + a, •) belongs to C^iR) and: 

^bxr^NSA/^ : „ M|2 



dx||e''^Gf*(x + a,-)lr < +oo. (71) 
We are now in a position to state and prove the main theorem of this section. 

Theorem lllli l: Let Va and Qa be defined, respectively, as in ([5T]) and (fH|) : let 6t and q solve 
the linear system ([57]) and 6*^ be the solution of Eq. ([52]) . Finally, let cpt = %4>-, with (j) £ >C^(M) 
and Tt defined as in (j6ip . Then the following three statements hold true with probability 1: 

1. Tt : £^(IR) — > £^(]R) defines a bounded operator for every t > (72) 

2. G ^b(IR) ^ 0f is a topological strong solution of Eq. ([6]) (73) 

3. G /:2(M) ^ lim||0i-(/)|| =0. (74) 

Proof of statement 1. Let 4> belong to £^(M); since also Gf^^(x,-) belongs to £^(M) for any 
X G M and t > 0, Holder's inequality implies that G^'°^{x^-)(f) belongs to £^(M); accordingly, the 
operator T^^'^ is well defined for any t > 0, and maps any £^(R)-function into a measurable 
function. By using Schwartz inequality together with relation (j67p . we have: 

j dx j dyGr\x,y)ct>{y) ' < Uf J dx WGf^ (x , ■)f < +oo; (75) 

thus T^NSA ^ belongs to £^(M) for any (j) in and for any t > 0. 

In a similar way, since also Gf^^(x + a, •) belongs to £^(M) for any a G C and because of ([7T|) . 
one proves that VoX^^^ (f) belongs to £^(M) for any G £'^(M), for any complex a and for any 
t > 0, i.e. that ViJ^a) contains TZ{T^^^). Using once more the same inequalities and (fTTI) . one 
shows also that Qb Va <j) belongs to C^{R) for any (p in C^{R), fr any a, 6 G C and t > 0. 

Remark: Actually a stronger statement is true, as can be readily seen from the Gaussian form of 
the Green's function Gt of the operator 7^: For positive t it maps 

£2(M) to Schwartz space o^Mj. 

We shall need this information in the proof of statement 3. 

Proof of statement 2. Let us consider the vector ipt := %^^^4>, with (p S >C^(M). By construc- 
tion, ipt solves Eq. (|59]) . once one proves that the integration 

•yNSA 



dyGr\x,y)^iy) (76) 

can be exchanged with the first and second partial derivatives with respect to x and with the 
first partial derivative with respect to t. We note that the function G^^^{x,y)4>{y) satisfies the 
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following two properties: i) The function y Gi^^{x,y)(f){y) is measurable and integrable on M 
for any t > and for any x G R; ii) The first and second partial derivatives with respect to x and 
the first partial derivatives with respect to t are exists for any t > 0, x G M and y € M and can be 
bounded uniformly with respect to t and x. Accordingly, one can apply e.g. theorem 12.13 pag. 
199 of [3] to conclude that the operations of integration and differentiation can be exchanged. 

Having proved that (ft solve Eq. (f59l) , a direct application of Ito calculus proves that (pt , defined 
as in ([6T|) . is a topological strong solution of Eq. ([6|). 



Proof of statement 3. Let cp = (pQ E C^(M) be given. Since (pt solves Eq. (I6|) in a strong sense, 
it also solves the SDE in a weak sense; hence, using e.g. Eq. (1.1) of {2^, one has: 

lhn{^\cPt) = {iflM G C-(M). (77) 

We extend ([77|1 to the general case of G £^(]R). Being dense in C'^{W), there exist a sequence 
{ifn G C^(M),n G N} which approximates any (p G £^(M). By triangle and Schwarz inequality we 
get 

\{ip\(t)t) - {ip\(t)o)\ < \{ipn\(l)t) - {^n\(l)o)\ + \\^-^n\\\\(l)t\\ + \\^-^n\\\\(t>o\\- (78) 

The first term on the right-hand-side can be made arbitrarily small because of ()77p : the second 
and third term can also be made arbitrarily small by choosing n sufficiently large, while \\(j)t\\ can 
be bounded as it converges to ||0o|| for t — > 0, due to Eq. Q. This proves that: 

lun{^\4>t) = G £\R). (79) 

Statement 3 for test functions (j) G (M) now follows directly from Eq. Q , Eq. (I79p and observing 
\\4't — '^olP = ll'/'tP + ll'^olP ~ 2((^o|0t)^- It remains to extend the strong continuity of % from the 
subspace C^(M) to C'^iR). For this observe that for cp G C^(M) {\\(pt\\^)t>o defines a stochastic 
process with continuous paths and by Holevo's result (cf. Eq. Q) it is a martingale. For given 
/ G £^(M) choose a sequence {(p^)neN C C^(M), which converges to / in £^(M). Doob's inequality 
for submartingales implies that for all n, m G N, T > and A > 

q( sup |||^rf-|l^rf|>A| < YlEQ[|||<^?^f (80) 

\0<t<T J A 

We now show that 

hm Eq[|||(^?;||2- ||^-||2|] =0. (81) 

n,m— +00 ^ L ij 
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The elementary inequality 

I II, II 2 II, „™ II 2 1 ^ IW.rJ^W I ll,„'Tt|| M|,„n ,„m,|| 

implies that 

M\Mf - \WTf\ < MiMW + \WT\\)M - vTW] 



< (EQ[(ii^rii + iiv^ni)'])MiEQ[iiv'r-^rf])^ 

The right hand side converges to as n, m — > 00. Therefore the sequence of stochastic processes 
(||(/7j^|P)t>o is a Cauchy sequence in the complete metric space (D,d) of adapted processes with 
right continuous paths having left limits, where the metric d is defined as (see page 56 - 57 in 
for background concerning this topology) 



d(x,y) = f;;iEQ Lnfi, sup Kx-y),]) 



n=l 

^n||2 



Therefore (Hv'j'll )t>o converges locally uniformly in probability to a stochastic process. This 
stochastic process again has to be continuous almost surely, since a subsequence of (||95"|P)t>o 
converges locally uniformly with probability one. Since lim„^oo llvrP = almost surely we 

know that [0, cxd) 9 t i-^ is continuous, in particular lim^^^o = ll/ll almost surely and 

defines by the lemma of Fatou a positive continuous supermartingale. Therefore it has a unique 
decomposition ||/t|p = Mt — At, where {Mt)t>Q is a continuous martingale and {At)t>o is increasing 
process. In fact as we shall show now, the increasing process is identically 0, i.e. ||/t||^>o is a posi- 
tive martingale for every / G C'^{M.). For that we observed in the Remark above that for positive 
e the function almost surely belongs to the Schwartz space and in particular to the domain of 
the generator. By Holevo's result cited above (||7f_e/e||)t>e is a continuous martingale. Therefore 
At = ^ for t > and hence it equals almost surely. In order to ensure strong convergence 
linif^oo Wft — /II = we need only show that weak convergence holds, i.e. \\m.t^oo{(l>\ft) = (01/) • 
Observing 

mft) - m\ < Mft) - {<t>wi)\ + m^i) - {<t>wn\ + \{<t>wn - m)\ 

it suffices to show that for some T > lim„_»oo supj<2^ K^l/t) ~ {4^Wt)\ = 0- But supj<2- K^l/i) ~ 
{(l^\Vt)\ ^ ll'/'ll sup«7^ ll/t— (^"||. Therefore we need only establish that lim„^oo supj<j. ||/f—99"|| = 0. 



20 



This is done by a similar argument as above, namely we show that for every e > 



hm Q sup||/t-<^rr >e = 0, 



because then there exists a subsequence which is almost surely convergent to 0. But as we showed 
above (||5't|P)t>o is a martingale for every 5 G / 
we can again apply Doob's inequality as before. 



above (||5't|P)t>o is a martingale for every g G Hence — ¥'"|P)t>o is a martingale and 



Remark 1. The Gaussian form of the Green's function (|12p is a consequence of the fact that 
Eq. dni) contains terms which are at most quadratic in q and p. T his also implies preserves shape 
of initially Gaussian wave functions; in fact, as shown e.g. in 



2g 



32 



331, 



3a] , a state 



t(x) = exp [-atix - xff + ikfx + , 



(82) 



is solution of Eq. ([6]) provided that the two real parameters x™ , A;™ and the two complex parameters 
(Tj, Q satisfy the following stochastic differential equations: 

2ih 



\- — {otY 
m 



ft , , ^ r 

dxf = —kfdt + —— 



dt, 

dit - 2^/Xxfdt 



dtJ 



VX^ \dCt - 2VXxf 



dq^ = ( X(xT? + —crl + -^]dt + ^/XxT dit - 2\/Axf dt 



m 



dd 



h 

2m 



dt + \/A-^ xf \d^t - 2VXxf dt 
erf L ' . 



(83) 
(84) 
(85) 
(86) 
(87) 



In particular, the solution of Eq. (f83l) is at = (A/u) coth(t;t + k), where k sets the initial condition. 
These results will be useful in the subsequent analysis. 



IV. REPRESENTATION OF THE SOLUTION IN TERMS OF EIGENSTATES OF THE 

NSA HARMONIC OSCILLATOR 

We now turn to the problem of analyzing the long time behavior of the solution of the (norm- 
preserving) non-linear Eq. ([1]). The representation of the solution (pt of Eq. ([6|) in terms of the 
Green's function (jl2p is not suitable for controlling the long time behavior; it turns out to be more 
convenient to express (pt in terms of the eigenstates of the NSA harmonic oscillator, resorting to the 
connection which we previously established between Eq. ([6]) and (I59p . In this way, as we shall see, 
the collapse process will be manifest: the coefficients of the superposition will decrease exponentially 
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in time, the damping being the faster, the higher the associated eigenstate. Accordingly — when 
normahzation is also taken into account — in the large time limit only the ground state survives, 
which has a Gaussian shape. 

We first recall a few basic features of the Hamiltonian of the NSA harmonic oscillator, 

K,2 



H 



— in\q 

2m 



(88) 



which has been studied in particular by Davies in a series of papers 
recent book 53|]. The eigenvalues of H are complex and equal to: 



51 



52l ] and reviewed in his 



l-i 



UJr, 



and the corresponding eigenvectors are: 

(t>^''\x) = ^/Ie-'^'''/^Hn{zx), z'^ = {l-i) 



Xm 



(89) 



(90) 



where Hn{x) is the normalized Hermite polynomial of degree n. Since the argument of Hn in ()90p 
is complex, these eigenstates are not orthogonal; it can be shown that they are linearly independent 
and form a complete set, however they do not form a basis. As such, they can not directly used 
to expand an initial state into a superposition of the eigenstates of H. This problem can be 
circumvented in the following way, also discussed by Davies. 

It is easy to see that the sequences and form a bi-orthonormal system; one then 

defines the (non-orthogonal) projection operators: 



Pn(t> = 

which satisfy the relations: 

p p — X p II p I 



and 



lim 



In WPr,. 



2c, 



(91) 



(92) 



where c is an appropriate constant 
sense that 



53 1 . As we see, although the states are normalized, in the 



(t)i^){x)(t>^'^\x)dx 



-'n,mi 



(93) 



the norm of the projection operators Pn grows exponentially as n ^ +oo. Finally, the following 



equality holds true 



53]: 



r^NSA ^ ^g-(i+iV„t/2p^ fort>4c/, 
n=0 



UJ. 



(94) 
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A remarkable property of the above representation of the solution of Eq. (j59p in terms of the 
eigenstates of the operator ()88p is that it holds not for any t > 0, as one would naively expect, but 
only for t > Ac/uj. The reason is that the norm of the projection operators Pn grows exponentially 
with n, so one has to wait for t to be large enough in order for the term g"""^*/^ to suppress 
the exponential growth of the projectors. From a physical point of view, recalling the discussion 
of sec HH since the constant c is of order 1 53|] and w ~ 5.01 x 10^^ sec~^, we see that the 
representation (j94p holds true only in part of the classical regime and in the diffusive regime, 
which is the one we are interested in studying now, but not in the physically more crucial collapse 
regime. 

We now apply the above results to our problem; we will first proceed in an informal way, and 
at the end we will prove the relevant theorems. Let (j) G /^^(M); then, according to ([6T]) and ([M|) : 



n=0 



(95) 
(96) 



n=0 



where a„ = {(^\(p^'^>*) (see Eq. (i9T]l ). while the two real parameters xt, kt and the complex parameter 
7t are defined as follows: 



xt = bf + bl- (2/mw)4 + {io/2VX)Ct, 
kt = {muj/h)bl + {l/h){cf -cl) + VX^t, 

2 



It 



-(1 - i){muj/Ah){bi - xi) + {i/h)et 



(97) 
(98) 
(99) 



By resorting to Eqs. ()57p and (j62l) . and after a rather long calculation, we obtain the following set 
of SDEs for these parameters: 



dxt = —ktdt + \ — 
m \ m 



dkt 



Xxi + - 



d^t — 2^/Xxtdt 
2^/Xxtdt , 

dt + VXxt d^t — 2VXxtdt 



h —2 iO 

k, H 

2m * 4 



dt 



VXxt d^t — 2\/~Xxtdt 



(100) 
(101) 
(102) 

(103) 



the initial conditions are: xq = /cq = To = 0. Note that these equations are equivalent to (I84l) - (j8/ 
with cjt = (Too = X/v = /2, Xt = xf, kt = k^ and 7^ = q + (1 + i)a;/4; as a matter of fact, the 
above equations describe the time evolution (according to Eq. ([6])) of the ground state of the NSA 
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harmonic oscillator, which is: 



(l)^{x) = exp 



-2^/ _ x2 --r l + i ' 
-y(x-xt) +iktx + jt 



(x) = (104) 



As we shall prove in the next section, this is the state to which — apart from normalization — any 
initial state converges to, in the long time limit, hence the name (pf. 

As we see, due to the stochastic part of the dynamics, the argument the Gaussian weighting 
factor and that of the Hermite polynomials of Eq. (|96p are different functions of time, while for 
analyzing the long time behavior of the wave function, it is more convenient that both arguments 
display the same time dependence. We thus modify the argument of the Hermite polynomials, to 
make it equal to that of the weighting factor. To this end, let us define Q = xt — h] we can then 
write: 

Hn[z{x-bt)] = ^ Hn[zix-xt) + zQ] 
/7r2"-n! 

1 \ - / n 



^i;(:)p^c.r-"^,„w.-.,) 



= E r-?; ' ,A ^^or-'"n,Mx-x,)], (105) 

— n V (n — m) 

where is the standard (not normalized) Hermite polynomial of degree m; in going from the first 
to the second line, we have used property (lA2p . Resorting to the above relation, we can rewrite 
Eq. ([96]) as follows: 

_ +0O 

^^(x) = e^^*^+^*-(i+'>*/^ a^'"'e-(i+^)"^^*/2^('")(x-Xf); (106) 

m=0 

the functions ^'^^^ are the eigenstates defined in (j90p . while the time dependent coefficients aj™' 
are defined as follows: 

where we have introduced the new quantity = e~(^"'"*^'^*/^Cj. 

Eqs. (llOOp and (|107p represent the two main formulas, which we will use in the next section to 
analyze the large time behavior. Before doing this, we need to set these formulas on a rigorous 
ground; we will do these with the following two lemmata. 

Lemma [iVll: Let </> G C?{^ and a„ = (0| with defined as in ([90]). Then the series (fTOT]) 

defining aj™' is a.s. convergent for any m and any t > 0. Moreover, one has the following bound 
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on the coefficients: 

+00 



V A;'^ 



fe=o 

wliere ^ is a constant independent of the Brownian motion ^, 



Proof: Because of ()92p . there exists a constant Ci such that: 

|a„| < ||<A||||<A^")|| = ll-^^^^ll < Cie"^ (109) 
Secondly, using Stirling formula, there exists a constant C2 such that: 

C2"iV2^n"e-" < n! < C2V2^n"e~", (110) 
for n > 1; we can then write the following estimate: 



m\k\ ~ V2vr V mk'^ ^m/2g-m/2^fcg-fe 

< ^2 ^-fc(lnfc-2)/2+m/2. /j^^^j^N 



in the second line, we have used the inequality (A; + m) ln(A; + m) < kink + mlnm + k + m. Using 
Eqs. (llOOp and (lllip . we have the following bound: 



m 



ak+m , (v2zCt 

Vmikl 



< (112) 



The cases /c = and m = can be treated separately, giving the same bound, with the only possible 
difference of an overall constant factor. This proves convergence of the series defined in (|107p and 
the bound (fTOSll . 

Theorem llVi l: Let the conditions of Lemma ITVl l be satisfied; let moreover 
where Ct = — bt with xt and bt solutions of Eq. (jlOOp and (j57p . respectively. Then the series 
defined in (jl06p is a.s. norm convergent for i > f = (4c + l)/w. In addition, the following equality 
holds true: 

+00 

Tt(f) = e^tx+^t~{\M)^tl^ ^ a^'"'e-(^+*)"'^*/V^™Ha; - xt), t > t, (113) 

m=0 

where % is the evolution operator associated to the Green's function ()12p . 
Proof: According to ([92]) and (jlOSp . one has: 

from which the conclusion follows. Comparing the two expressions of Eq. ()6ip and Eq. ()106p when 
the initial state i;^ is an eigenstate (p^""^ , we see that they coincide on the dense subspace of all finite 
linear combinations of (i^"^ , and hence on the whole of (M) . 
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V. THE LONG TIME BEHAVIOR 

We are now in a position to study the long time behavior of the solution of Eq. ([T]) . Looking at 
expressions (I106P for the solution (j)t and (I107P for the coefficients it should be clear what the 
long time behavior of the normalized solution ipt = <Pt/\\(pt\\ is: whatever the initial condition, at 
any time t > the wave function (pt picks up a component on the ground state (p^^^x — xt), since 
af^ 7^ as long as at least one of the coefficients is not null, which is always the case. Eq. (jl06p 
on the other hand shows that each term of the superposition has an exponential damping factor, 
which is the bigger, the higher the eigenvalue. Accordingly, after normalization, only the eigenstate 
with the weakest damping factor survives, which is the ground state. Hence we expect that the 
general solution of Eq. ([1]) converges a.s., in the large time limit, to the ground state (p^'^^x — xt), 
which is a Gaussian state. That this is true is proven in the following theorem. 

Theorem [Vll: Let (pt be a strong solution of Eq. ([6]) that admits, for t > t a representation as 
in ()113p . Let ^pt = <Pt/\\4't\\ (when \\(pt\\ 7^ 0), which can be written as follows: 

+ 00 — (m) 



i{ktx+^l-ujt/4:) 



(115) 



with: 



_(0) 

J(fctx+7j-a;t/4) 



(po{x - Xt), 



+ 00 



m=0 



(116) 
(117) 



Then, with P-probability 1: 



lim 

t—*oo 



0. 



(118) 



Note that, apart from global factors, tp^ is the ground state of the NSA harmonic oscillator, 
randomly displaced both in position space as well as in momentum space. 



Proof. According to Eq. (jllSp . all we need to prove is that, with P-probability 1: 



lim 



+ 00 (m) 

^1 



Resorting to (jll4p . one can write the following bound: 



(119) 



+ 00 (m) 

V^e-(i+^)™-*/V™(x-x,; 

m=l 



< Ci — 



n 1- e-^(*-*) 



(120) 
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thus all we need to set is the long time behavior of rt and Nf. Lemmas IVl 1 and|Vl2 (see Eqs. (jl2ip 
and (I126p ) state that, with P-probability 1, rt converges asymptotically to a finite and non-null 
random variable, while Nt converges to a finite random variable. From these properties, the 
conclusion of the theorem follows immediately. 

In the remaining of the section, we prove the required lemmas. 

Lemma [Vll: Let n be defined as in (jllTh . Then, with P-probability 1, 



lim rt = Too finite and not null. 

t—rOO 



(121) 



Proof. According to Eq. (|106p and (I117p . the following equality holds: 



(122) 



resorting to the stochastic differentials Q and (|102p for 



and 7^ respectively, one can write 



down the following stochastic differential equation for r|: 



dr^ 



2^{{q)t - xt) dit + 4A {x1 - {q)t Xt) dt 



1 1 



1. 



(123) 



By using relation (jlip . the above equation can be re- written in terms of the Wiener process Wt as 
follows: 



dr^ = {{q)t - Xt) dWt + 4A {{q)t - Xtf dt 



^0 — J-) 



(124) 



whose solution is: 



exp 



2VA / {{q)s-Xs)dWs + 2\ / {{q)s-x,fds 
Jo JO 



(125) 



The crucial point is to establish the behavior of the difference {q)t — xt between the mean 
position of the general solution ipt and the mean position of the "asymptotic" state ip^. Since 
ipt converges to ^pt°, we expect {q)t — xt to vanishes asymptotically. That this is actually true 
with P-probability 1 is proven in Lemma W\3 (see Eq. (I129|) ). where indeed it is shown that the 
convergence is exponentially fast. This fact, together with (I125p . concludes the proof of the lemma. 

Lemma [Vl2: Let Nt be defined as in (fTOSi) . Then, with P-probability 1, 



lim Nt = Nao finite. 

t— »oo 



(126) 
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Proof. Looking back at Eq. (jlOSp . we see that in order to prove this lemma it is sufficient to 
show that Ci tends to a finite hmit as t — > oo, with P-probabihty 1. According to our previous 
definition, C,^ is equal to: 

It = e-^^+'^-'l\xt-htY (127) 

Eqs. (f57l) and (f97|) . together with the change of measure ([TT]) . lead to the following stochastic 
differential equation for in terms of the Wiener process Wt'- 

dWt + 2VX{{q)t - xt)dt^ , Co = 0. (128) 

Once again, the large time behavior of {q)t — xt (see Eq. (jl29p ) yields the conclusion of the lemma. 

Lemma [Vl3: Let {q)t = (V'tkl^t) and xt defined in ([97]) . Then, with P-probability 1: 

ht ^ {q)t-xt = 0(e--*/2)_ (^29) 



dCt 



^ _^-(l+»)a;f/2 



2\/A 



Proof. Let us consider the Gaussian solution of Eq. 

0f(x) = Gt(x,0) = i^texp 



at 2 — - ' 
—X + atx + ct 



at , _c X 2 -tG 

- xX) + iktx + Ct 



where Gt{x,y) is the Green's function defined in ()12p and 



a 



a; 



kt — «i 5 



a; 



ct = ct + ^{xf?. 



(130) 
(131) 

(132) 



Note that xf is the mean position of the Gaussian state cpf, while hk^ is its average momentum. 
Obviously we can write: 



ht = {{q)t-xf) + {xf -Xt); 



(133) 



lemma |BJ1 proves that {q)t — xf has the required asymptotic behavior (see Eq. ()Bip ). so all we 
need to show is that also xf —Xt behaves as required. Lemma [Bjl was first proven in [33]; for 
completeness, we reproduce it in Appendix[Bl adapting it to our notation. The proof of the lemma 
is instructive because it makes clear why it is convenient to analyze {q)t — xf separately from 
xf - Xt. 

By letting the ground state of the NSA harmonic oscillator evolve according to the Green's func- 
tion Gt{x,y), one can express xt in terms of the functions (|13p ~ (ll8p : a straightforward calculation 
leads to the following result: 



(p--iW-(Aj 

\at + aoo/ 



(134) 
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where aoo = linif^oo a* = 2A/f. By inspecting expressions ([Ml) and (fT5]) . we recognize that 
p^^ — 1 = 0(e~'^*) and \f3t\ = 0(e~'^*/^), thus in order to prove the lemma aU we have to do is to 
control the long time behavior of a^, which in turn sets the asymptotic behavior of bt through (fTT]) . 
Inverting Eq. (jl32p we get: 



at 



-G , -T^ 

at Xf + ik^ , 



(135) 



thus we can control at by controlling and ki . These two quantities, being the average position 
and (modulo h) average momentum of the Gaussian solution (jlSip , satisfy the stochastic differential 
equations (j84|) and (|85|) . with at/ 2 in place of ctj. By using the change of measure (|lip . we can 
re-express these equations in terms of the Wiener process Wt as follows: 



dx? 



dk 



h-G 
— k 

m 



a 



dt + ^ dWt, 



2^^ftdt - Vx^dWt, 
a+ ay- 



(136) 
(137) 



with ft = {q)t — xf. By integrating the second equation, by using the strong law of large numbers 
applied to Wt, Eq. (|Bip for ft and the fact that at has an asymptotic finite limit, one can show 
that, with P-probability 1, the process k^ grows slower than t^, for t — > oo. By integrating now 
the first equation, and by using the same properties as before, one can show that xf grows slower 
than t^, for t — > cxd and again with P-probability 1. According to Eq. ()135p and ()17|) . we then have, 
with P-probability 1: 



at = o{t'^) as t — > oo. 



lim bt 



bco finite. 



(138) 



This proves that xt — xf has the required asymptotic behavior, hence the conclusion of the lemma. 

In this way we have proven that any initial state is P-a.s. norm convergent to the Gaussian 
state ()116p . which can be written as follows: 

v2 



vr 



- — {x-xt) + iktx + i - jt 



which has a fixed finite spread both in position and in momentum, given by 

A, 



V mio 



hmuj 



(139) 

(140) 
(141) 



This corresponds almost to the minimum allowed by Heisenberg's uncertainty relations, as A^Ap = 
hl\pi. Note also that, the more massive the particle, the smaller the spread in position of the 
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asymptotic Gaussian state: this is a well known effect of the localizing property of Eq. ([T]). Finally, 
Eqs. (jlOOp and (llOip . together with the change of measure (jlip . tell how the average position Xt 
and momentum hkt evolve in time, as a function of the Wiener process Wt- 

UJ 



dxt 
dkt 



— ktdt + ijjhtdt + — i=dWt, 
m 2VA 



2Xhtdt + VXdWt, 

which imply that there exist two random variables X and K such that 

Xt = X + —Kt + VX— [ Wsds + \ —Wt + 0(e 
m 1^ Jo \ m 

h = K + VXWt + Oie"^'^^). 



-ut/2\ 



(142) 
(143) 

(144) 
(145) 



These parameters fully describe the time evolution of the Gaussian state (|139p . 



VI. CONCLUSIONS AND OUTLOOK 

In section |II] we have spotted three interesting time regimes during which the wave function, 
depending on the values of the parameters A and m, evolves in a different way. In the central 
sections of this paper we have analyzed the long time behavior, which pertains to the third regime, 
the diffusive one. There are many other properties of the solutions of Eq. ([1]) which deserve to be 
analyzed, and in this conclusive section we would like to point out a number of interesting open 
problems. 

I: Collapse regime. Let ^ be the length which discriminates between a localized and a non- 
localized wave function, i.e. such that, defining with /X^q the spread in position of a wave function 
■0, we say that V is localized in space whenever /X^q < £. In our case, we must take £ > y^T™^) 
where sjfilrrujj is the asymptotic spread (see Eq. (I140j) ). 

Problem I.l: collapse time. Let be the solution of Eq. ([T]), for a given initial condition 
'(/' G £^(IR) such that X!^q > i. Let us define the collapse time Tj^ol ^ the first time at which the 
wave function is localized in space: 

r^oL := min{t : X'^'q < £}. (146) 

Question 1. 1.1: How is TcoL distrib ute d, as a random variable? In particular, is it finite with 
P-probability 1, as we expect it to be [so!]? What are its mean Ep[T|^ql] and variance Vp[T^ql]? 
Question 1.1.2: How does T^ql depend on the initial spread X^q, as well as on the parameters 
A and m? 
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Question 1.1.3: What is the probabihty that, for t > Tqq^, the wave function de-locahzes in space, 
i.e. acquires a spread greater than £, namely A^q > £ + e, where e is an arbitrary positive quantity? 
This kind of analysis is important because it gives a measure of how stable the localization process 
is. 

Problem 1.2: collapse probability. Let ip := i^t, for t = Ep[Tcol]- Let x := {^p\q\'^l^) be the 
position of the wave function at the average time at which it is localized in space. 
Question 1.2.1: How is x distributed as a random variable? 

Question 1.2.2: Let p{x) be the probability density of re; let be the collapses probability 

density given by the Born probability rule. When does it happen that 

d{x) := \p{x) - \i:{x)f\ < 6, (147) 

where 6 is an appropriately small number? How does this depend on the values of the parameters 
A and m? 

II: Classical regime, in the classical regime, the wave function is expected to move, on the 
average, like a classical free particle. 

Problem II. 1: classical motion. Let and be the (quantum) average position and mo- 
mentum of tpf Let t > Tf^oL- 

Question II. 1.1: How are q^ and p^ distributed, as random variables? In particular, what are 

their mean Ep[qJ, ]Ep[gJ and variances Vp[gJ, Vp[gJ? 

Question II. 1.2: How do they depend on the values of A and m? 

Question II. 1.3: Let T^jp be te time at which the motion departs from the classical one 

T^ip := min{t : Vp[gJ > Ag V Vp[pJ > Ap}, (148) 

where Aq and Ap are suitable parameters measuring the fluctuations of the position and momentum, 
respectively, of the wave function. How does T^ip depend on the parameters of the model? 

Ill: Diffusive regime. This regime begins after T^jp, and it has been analyzed in this paper: as 
we have seen, the wave keeps diffusing in the Hilbert space, eventually taking a Gaussian shape, 
as described in Sec. FVl 
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APPENDIX A: PROPERTIES OF HERMITE POLYNOMIALS 



We list here the main properties of Hermite polynomials, which are used in the paper. The 
primary definition of the Hermite polynomials is 



Ln/2J 



Hn{z) = n! ^ 



m ^\n—2m 



m=Q 



m!(n — 2m)! ' 



where z is any complex number. These polynomials satisfy the following addition rule 

m=0 ^ ^ 



(Al) 



(A2) 



When the argument is real {z = x G M), they form an orthogonal set with respect to the weight 
exp[— x'^]; the normalized Hermite polynomials are: 

1 



Hn{x) 



Nr. 



Hn{x), Nn = V V^2"n!. 



(A3) 



APPENDIX B: LEMMA 



Lemma Ell: Let <p e £^(M), \\(t>\\ = 1 and let (pt = Then, with P-probability 1: 

ft = (g)t-xf = 0(e-*/2), 
where {q)t = {'4't\<l\i^t) , and has been defined in (1132p . 



(Bl) 



Proof. Using the expression (I12p for Gt{x,y) together with Schwartz inequality, we can derive 
the following bound on cpt- 



< \Kt\ 




exp 



,^Pt-H. 



* x2 + 2 I a? + 



X + 2cf + 



(B2) 



^ Pt \^ Pt J ~ 2\ pt ^ 

which holds for any t > 0. The above inequality implies that it is sufficient to consider </> G £^(M) 
such that: 



{x)\ < Ce 



(B3) 



where C and A are random variables. A direct calculation leads to the following expression for the 
quantum average {(j)t\q\4>t)'- 



h\q\(pt) 



■ exp 



■ exp 




exp 



dyidy2<t){yi)(p{y2 



Ptyi + Pty2 ^ at 

2af 



2 V"* 2af 



Of 



2af 



2/2 



(B4) 
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As we shall soon see, all exponential terms in the above expression can be controlled. The 
crucial factors are the two within brackets: the first term decays exponentially in time, since 
(it = 0(e~'^*/^), while at has a finite asymptotic limit; the term af /af , instead, does not decay in 



time (see the discussion in connection with the proof of lemma IVl3). Since 



is equal to the 



expression (IB4p without the terms in square brackets, and because of (I132p . we have that 



(B5) 



\Kt\ 



■ exp 



• exp 




a exp 



2cf + 



a; 



at 



2af ) ^1 2 



at 



bt + ^]yi + {bt + '^' 



dyidy2<t){yi)<P{y2) 



2af 



2af 



yl 



y2 + TTTir ^1^2 



a 



2aJ 



(B6) 



According to the discussion above, we expect the quantity ft 



to decay exponentially in time, 



as we shall now prove; this is the reason why, in proving lemma |Vl3, it was convenient to split the 
difference ht as done in Eq. (I133p . 

Using the inequality yiy2 < {yl + yi)/^ we can write: 



t\\m\ 



af 



< i^\Kt\ J-^exp 




2cf + 



a. 



(B7) 



dyidy2 |(/>(yi)l 10(^2)1 (|yi| + |y2|)5'(yi)5(2/2),(B8) 



with: 



9{y) 



exp 



1 



at 



(A 



R\2 



a; 



y' + [bU 



a? 



(B9) 



Next, by using the inequality g{yi) +5(2/2) < {oiyi)^ + and the symmetry between yi 

and y2, we have: 



l/t 



< 



2af 



\Kt\ 




exp 



■R^2 



af 



dyidy2 |</'(2/i)ll0(y2)|(|2/i| + \y2\)giyi 



Now, a direct computation shows that 



\Gt{;: 



dx\Gt{x,': 



\Kt 




-exp 



2cf + 



a\ 



9{y?; 



(BIO) 



(Bll) 



the key point is that, since Gt{x,y) solves Eq. ([6]), then 



is a positive martingale with 



respect to the measure Q, for any value of y; we call Mar(Q(t, y) this martingale. We can then write: 



\ft 



< 



m 

2af 



dyidy2 |</'(yi)||0(2/2)|(|yi| + \y2\)MarQ{t,y) 



-2^/ dye-^y\Ai\y\+A2)MarQ{t,y), 



(B12) 
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where Ai and A2 are suitable constants. In going from the first to the second hne, we have 
used (IB3p . The quantity 

^jdy e-^y" {Ai\y\ + ^2)MarQ(i, y) (B13) 

is another positive martingale with respect to Q, which we call MarQ(t). We arrive in this way at 
the inequality: 

\ft\ < \(3t\ (B14) 

Since MarQ(t) is a positive martingale with respect to Q, then Marp(t) = MarQ(i)/||i;^>t|p is a 
positive martingale with respect to P which, by Doob's convergence theorem, has a P-a.s. finite 
limit for t — > +00. The conclusion of the lemma then follows from Eq. (jlSp . according to which 
= 0(e-^*/2). 
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